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How fast do Jupiters grow? Signatures of the snowline and 
growth rate in the distribution of gas giant planets 
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ABSTRACT 

We present here observational evidence that the snowline plays a significant role in 
the formation and evolution of gas giant planets. When considering the population of 
observed exoplanets, we find a boundary in mass-semimajor axis space that suggests 
planets are preferentially found beyond the snowline prior to undergoing gap-opening 
inward migration and associated gas accretion. This is consistent with theoretical 
models suggesting that sudden changes in opacity - as would occur at the snowline 
- can influence core migration. Furthermore, population synthesis modelling suggests 
that this boundary implies that gas giant planets accrete ~ 70 % of the inward flowing 
gas, allowing ~ 30 % through to the inner disc. This is qualitatively consistent with 
observations of transition discs suggesting the presence of inner holes, despite there 
being ongoing gas accretion. 

Key words: planets and satellites: formation — Solar system: formation — stars: 
pre-main-sequence — planetary systems — planetary systems: formation — planetary 
systems: protoplanetary discs 



1 INTRODUCTION 

Since the discovery, in 1995, of the first extrasolar planet 
around a main-sequence star (|Mavor fc Queloj|l995l ) , a fur- 
ther 759 such planets have been detected. Many of these 
planets have masses similar to that of Jupiter and hence are 
generally regarded as gas giants. The standard model for 
the formation of th ese planets is the core accretion model 
( Poll ack et al.|[l99c3 ) . In this model, micron-sized dust grains 
grow to form kilometre-sized planetesimals that then coagu- 
late to form planetary mass bodies that, if sufficiently mas- 
sive, may gravitationally attract a gaseous envelope if gas is 

still present in the disc. 

Populatio n synthesis models (|lda fc Lin! 1 20041: 

lAlibert et al.l [2005) have generally been successful in 
reproducing the properties of the observed exoplanet 
population. However, while these models have illustrated 
how the overall process leads to a population consistent 
with that observed, they have not specifically quantified 

any individual parts of the process. 

It has, however, recently been noted l|Wright et al.l 

2009) that the distribution in semimajor axis of the ob- 
served exoplanets shows a peak at ~ 1 AU. It has been sug- 
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gested that this peak may be a consequence of disc disper - 
sal through photoevaporation llAlexander fc Pascuccill2 012). 
Some population synthesis models ( Mordasini et al.l I2012T ) 
have also managed to reproduce this peak in the semima- 
jor axis distribution and suggest that it is an imprint of the 
snowline. The ices that collect on small dust grains subli- 
mate inside the snowline, where the temperature in the disc 
exceeds 17 K. The snow line radius is 2.7 AU for a solar- 
mass star (Hayashi 1981) and scales with M 2 . 

During the gas giant planet formation process, the 
rocky /icy core is expected to go through a p hase of rapi d 
inward migration, known as type I migration (IWardlll997l ). 
Quite how these cores survive to form gas giant planets is 
still uncertain. This migration may, however, be strongly 
affected by regions with a su dden change in opacity, a s 
would occur at the snowline l|Menou fc Goodman! 120041 ) . 
To investigate this we consider the distribution of planets 
in semi-major axis space. In the left-hand panel of Fig- 
ure [1] we plot planet mass against semimajor axis, while 
in the right-hand panel we plot planet mass against semi- 
major axis with the semimajor axis normalised with re- 
spect to the predicted snowline of the host star. We con- 
sider 463 exoplanets that were first detected by radial ve- 
locity measurements. In the left-hand panel there appears 
to be an e xcess of plane t s beyo nd 1 AU th at has been high- 
lighted bv lWright et al.l {2009) and others (|Mordasini et all 
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Figu re 1. The left-hand panel shows planet mass against semimajor axis and illustrates an excess of planets beyond 1 AU l|Wright et al.l 
2009). The right-hand panels shows planet mass against semimajor axis with the semimajor axis normalised with respect to the snowline 
of each exoplanets host star. In the right-hand panel, there appears to be quite a well-defined diagonal boundary, across which there is 
a step increase in the density of planets. 
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Figure 2. Planet mass against semimajor axis with the semi- 
major axis normalised with respect to the snowline of each exo- 
planets host star. There appears to be a diagonal boundary, in 
mass-semimajor axis space, across which there is a step increase 
in the density of planets. The solid line shows the boundary that 
produces the the largest density contrast while the dashed and 
dotted lines shows two other solutions that lie on the extremes of 
the bootstrapped 49% confidence contour in 8 — r\ space. 



I2012I : lAlexander fc Pascuccill2012l ) . In the right-hand panel, 
however, there appears to be, for planet masses below a 
few Jupiter masses, quite a well-defined diagonal boundary 
across which there is a step increase in the density of planets. 

To characterise the diagonal boundary seen in the right- 
hand panel of Figure Q] we assume that it has the form 
Msini/Mj U p = S (a/a anov ) v . We consider a box with 0.2 ^ 
M/Mj up < 7 and 0.07 ^ a/a snow < 1 and determine the 
values of <5 and r\ that produce the biggest density contrast. 
The solid line shows the boundary that produces the largest 
density contrast while the dashed and dotted lines shows two 
other solutions that lie on the extremes of the bootstrapped 
49% confidence contour in S — r\ space, discussed in more 
detail in Section [Al 

If a planet is sufficiently massive, it can open a gap 



and migrate inwards via what i s known as type II migr ation 
|Goldreich fc Tremainei 1 198(1 iLin fc Papaloizoul Il986h . It 
can also continue gro wing by accreting some of the gas flow - 
ing though the gap l|Lubow. Seibert fc Artvmowicj[l999l ). 
During this phase, planet growth typically follow s a diagonal 
line in log M - log o space l|Mordasini et al,ll2012l ). We, there- 
fore, assume that the diagonal boundary in the right-hand 
panel of Figure [TJ and characterised in Figure [2] is a con- 
sequence of this migration and gas accretion. The diagonal 
lines in Figure[2]have the form Msini/Af.j up = S (a/a S now)' J 
and so if a = a snow , M = 5 Mj up , which we interpret as 
suggesting that a planet forming at the snowline reaches a 
mass of 8 Mj up before starting to migrate inwards via type II 
migration and growing via associated gas accretion. For the 
best fit line (solid line), this corresponds to M a = 0.01 M.j up , 
while for the dashed line it is M — 0.078 Mj up . These are 
both consistent with the mass at which we would expect type 
II mi gration to start operating (|D'Angelo. Henning fc Klevl 
120031 ). The dotted line gives M a < 0.01 M Jup which may be 
unphysically small. We can, however, then use this to pre- 
dict where each planet was prior to the start of this process. 
If this process typically starts when the planets reach a mass 
of M , then from its current mass, M p , and semimajor axis, 
a, its initial semimajor axis would be a a = a(Af /M p ) 1// ' 7 . 
We thus determine the initial distribution of these planets, 
as shown in Figure [3] where the initial semimajor axis is 
normalised to the snowline and the line styles correspond 
to those in Figure [2] Figure [3] shows a definite jump at the 
snowline, with the largest density contrast occuring, as ex- 
pected, for our best fit line. 

The above suggests that prior to the final stages of 
planet formation (inward type II migration and final stage 
gas accretion), planets are preferentially found beyond the 
snowline of their host star. In this paper we, therefore, com- 
bine a realistic model of the evolution of protostellar discs 
with models of the growth and migration of gas giant plan- 
ets to establish if this intepretation of the diagonal feature 
in Figure [2] is consistent with theoretical models of planet 
migration and growth. 
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Figure 3. Distribution of the initial semimajor axis, a , of the 
exoplanet population with a normalised with respect to the ex- 
pected snowline of the host star. In this case a a is determined 
by assuming that growth follows the diagonal line in Figure [2] 
and hence that a a = a(M /A/ p )~ 1 < / ' ! with M a = M Jup . The 
line styles correspond to those in Figure [2] and the largest density 
contrast occurs, as expected, for the best fit line in Figure [2] 



2 BASIC MODEL 
2.1 Disc model 

We assume that the disc is axis ymmetric and use the 
standard one-dimen sional equations |Lvnden-Bell fe Pringld 
Il974l ; iPringle] Il98lh to evolve the surface density, £(r,t). 
The surface density evolution is largely determined by the 
kinematic viscosity, v, and by mass loss through a disc 
wind. We a ssume that the viscosity has the form of an 
a viscosity (jShakura fc Sunvae v 1973) so that v = ac a H, 
where c 3 is the disc sound speed, H is the disc scale- 
height, and q << 1 is a parameter that determines the 
efficiency of an gular momentum transport. It has recently 
been suggested (|Owen et al.ll2010r ) that x-rays are the dom- 
inant driver of photoevaporation, and so here we imple- 
m ent the x-ray photoionizat i on m odel described in detail 
in lOwen. Ercolano fc Clarke! |201ll ). 

We ran 100 disc models and selected the central star 
mass randomly between M* = 0.8 Mq and M* = 1.2 Mq. 
We don't explicitly include the planets in these disc mod- 
els, but instead use the range of disc models to later evolve 
planets with a large range of different initial conditions. We 
assume the disc extends from r = 0.1 AU to r — 50 AU with 
an initial surface density profile of E oc r~ . In each simu- 
lation the initial disc mass is 0.25M*. The initial disc mass 
is therefore quite high and s uch discs are likely to be self- 
gravitating. Within 50 AU (Rafikov 2005), it is expected 
that such discs will achieve a state of quasi-steady ther- 
mal equilibrium with dissipation due to the gravita tional 
instability balan ced by radiative cooling |Gammiell200ll). It 
has b een shown (|Balbus fc Papaloizoull 19991 ; Lodato fc Rice! 
|2004| ) that the gravitational instability then acts to trans- 
port angular momentum in a man ner analogous to viscou s 
transport. As described in detail in Rice fc Armitaga (2009) 
(see also IClarkei 120091 ; IZhu. Hartmann fc Gammiei 120091) . 
this can then be used to determine the effective value of 
a and, hence, the kinematic viscosity, v. 

If, however, the effective gravitational a is less than 



0.005, we assume that another transport mechanism, such as 
the m agnetorotational instability (MRI) |Balbus fc Hawlevl 
Il99ll) will then dominate and we set a = 0.005. In addi- 
tion, we also assume that irradiation from the central star 
sets a radially dependent minimum temperature in the disc 
(|Havashilll98ll ). Although all of our disc models start with 
the same disc-to-star mass ratio, the variation in x-ray lumi- 
nosity (from 5 x 10 2S erg s _1 to 10 31 erg s _1 ) produces a wide 
range of d ifferent disc lifetimes, largel y consistent with that 
observed (jHaisch. Lada fc Ladall200ll ). We therefore have a 
set of disc models that can self-consistently evolve the sur- 
face density from the early stages, when the gravitational 
instability is likely to dominate, through to the later stages 
when an alternative transport mechanism, such as MRI, will 
dominate and that also includes the late-stage dispersal due 
to photoevaporative mass-loss. 



2.2 Core growth and migration 

Planet formation is thought to occur through the initial 
growth of micron-sized dust grains to, ultimately, kilometre- 
size planetesimals. These planetesimals then continue to 
grow with, typically, the largest in any region dominat- 
ing and undergoing what is known as oligarchic growth 
(|Kokubo fc Idall998h . If these oligarchs can grow sufficiently 
massive (~ 10 M®) then, if there is still gas present in the 
disc, they can gravitationally attract a gaseous envelope and 
can rapidly grow to become a gas giant planet. 

Once a planetary mass body has formed, it can ex- 
change angular momentum with the surrounding disc ma - 
terial and migrate radially (|Goldreich fc Tremaind fl980l ). 
There a number of different migration scenarios. Low-mass 
planets (M p i ~ 1 Me) generate a linear disc response and 
migr ate inwards through what is known as type I migra- 
tion (Ward 1997). High-mass planets are thought to open a 
gap in the disc gas and migrate through type II migration 
(|Lin fc Papaloizou|[l9 86). A third type of migration, known 
as ty pe III migration, may occur for intermediate mass plan- 
ets (|Masset fc Papaloizoul 120031 ). In the type III regime, 
corotation torques gen erate very rapid migration w hich can 
be inward or outward (|Masset fc Papaloizoul [20031 ). 

The migration rate for planets migrating in the 
type I regime increases with planet mass and an- 
alytic jTanaka. Takeuchi fc Ward] |2002| ) and numerical 
ijBate et al l 120031 ) estimates suggest that the timescale 
should be less than typical disc lifetimes. These plan- 
ets should therefore migrate into the central star prior 
to becoming massive enough to enter the slower type 
II migration regime. Recent numerical simulations sug- 
gest that the type I migration rate can, however, 
be significantly reduced if the disc thermodynamics is 
treated more realistically JPaardekooper fc Mellerrial 120081 ; 



iPaardekooper. Baruteau fc Klevl 201ll ). Recent population 
synthesis models therefore typically assume that the an- 
alytic ty pe I rate is reduc ed by at least a factor of 
about 30 (|Alibert et al.ll2005l ). These models are now quite 
successful at reproducing the observed characteristics of 
the exoplanet population dMordasini. Alibert fc Benzj|2009l ; 



lAlibert. Mordasini fc Benzll2011 



Furthermore, 



has 



been suggested 
(IMenou fc Goodman] |2004| ) that type I migration can 
be strongly influenced by sudden changes in the disc opac- 
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ity, such as may occur near the snowline. One of the goals 
of this paper is to investigate whether or not the possible 
infuence of the snowline on core migration is consistent 
with the observed distribution of extraolar planets. After 
type I and, potentially, type III migration have distributed 
the planetary cores throughout the disc, they then undergo 
continued growth through gas accretion and migrate via 
gap opening type II migration. We also want to compare 
the current distribution of exoplanets with a population 
synthesis-like model to determine if we can quantify type II 
migration and the growth that accompanies this migration 
process. 

2.3 Type II migration and gas accretion 

In our modelling we don't explicitly model the growth 
of the cores and the initial runaway gas accretion phase 
l|lkoma. Nakazawa fc Emorill200ol ; iBrvden et al.ll200d ) that 
occurs when the planetary core reaches the critical core mass 
lllkoma. Nakazawa fc Emoril |2000| ; iPapaloizou fc Terqueml 
Il999l ). The runaway gas accretion phase t erminates whe n 
the planet reaches the gas isolation mass (|Lissaueij[l987l ). 
We assume that the gas feeding zone is 2 Hill radii wide 
with the Hill radius of a planet with mass M p located at a 
in the disc given by 



th = a 



( M p 
V 3M» 



1/3 



(1) 



One can then show that the isolation mass is 

lso ~ (3M,)i/2 ' W 

As discussed in section |2~T1 we have a series of disc mod- 
els with central star masses between M 4 = 0.8 Mq and 
M* = 1.2 Mq and with a range of x-ray luminosities that re- 
sults in a range of disc lifetimes that largely matches that ob- 
served (|Haisch. Lada fc Ladall2~00ll ). The current exoplanet 
sample has more stars with mass M, > 1 Mq than with 
masses M 4 < 1 Mq and so we use the disc models with cen- 
tral star masses > 1 Mq twice as often as those with central 
star masses < 1 M . 

For each planet formation simulation, we assume that 
the planetary core forms and undergoes runaway gas accre- 
tion at a randomly chosen time between t = 1 Myr and 
t = 4 Myr. The disc model then gives us the mass of the 
central star and surface density and so Equation ((2]) can 
be used to determine the planet's isolation mass. The disc 
model also gives the time dependence of the disc viscosity, 
v. If sufficiently massive, a planet will open a gap in the 
disc dLin fc Papaloizoulll986l ). For a planet located at a, the 
width of the gap, A, satisfies (|Sver fc C larke 1995|) 



Qa 2 



(3) 



where q — M p /M t . For a gap to open, the disc scaleheight, 
H, must be less than the gap width A. A second criterion 
is that the gap width needs to be greater than the Roche 
radius, Rl, of the planet and hence 



A > Rl = q 



1/3, 



(4) 



Type II migrati on effectively ha s two regimes, the "disc 
dominated" reg ime ( Armitagc 2007) and the "planet dom- 
inated" regime (|Trilling et al.lll998l ). "Disc dominated" mi- 
gration occurs when the local disc mass is large compared 
to the mass of the planet. In this case the planet is coupled 
to the viscous evolution of the disc and the migration rate 
is independent of the mass of the planet. The radial velocity 
of the planet is then 

«W = ~. (5) 

If, however, M p > 2Ea 2 , the inertia of the planet reduces its 
radial velocity to 

3u Ea 2 



(6) 



If Equation ((4]) is satisfied and H < A then the planet is 
able to open a gap and migrate through type II migration. 



While migrating inwards, the planet is also able to accrete 
mass from the disc. For planets with mass of a few J upiter 
masses, it has been sugg ested (|Klev fc Dirksenl fcooe) that 
the planet could accrete at a rate comparable to gas ac- 
cretion rate through the disc (M 4 = 37r^E). We therefore 
assume that in the "planet dominated" migration regime 
the accretion rate onto the planet is 

M p = /33tti/E, (7) 

with B typically close to unity (|Mordasini. Alibert fc Benzl 
(2009) assume B = 1), and v and E coming from 
the self-consistent disc models. In the "disc dominated" 
regime we assume that this is reduced by a factor 
1/B where B = (2nT,a 2 ) /M p is essentially the ra- 
tio of the local disc mass to the planet mass. The 
reduction is motivated by the ma ss accretion effi- 
ciency for planet growth determined by IVeras fc Armitagel 
(2004) based on t he results of two-dimensional nu- 
merical simulations fLubow. Seibert fc Artvmowicz! 1 19991 : 
iD'Angelo. Henning fc Klevll2002r i. 

We can combine Equations [6] and [7] to show that, in 
the "planet dominated" regime log M p /M = log (a/a )~ 137 ' . 
The dotted line in Figure [5] has a slope steeper than — n 
which, given that B ^ 1, would appear to be unphysical 
if our simple interpretation is correct. We almost certainly 
need more data to determine, more accurately, the proper- 
ties of this boundary. 



2.4 Putting it all together 

We randomly select one of our completed disc models and 
start modelling the evolution of the planet at the stage at 
which the planet semimajor axis has been set by type I 
(or type III) migration and the planet has just undergone 
runaway gas accretion. The initial planet mass is given by 
Equation ([2]) with E taken from our chosen time-dependent 
disc model. Given the planet mass and semimajor axis, we 
use Equations (J3]) and Q to determine if the planet opens 
a gap in the disc. If not, we assume that it grows slowly 
(M p = 0.1M* with M, determined from our chosen disc 
model) and remains at its initial semimajor axis until it ei- 
ther does satisfy the gap opening criteria, or the gas disc 
has dispersed. If the gap opening criteria are satisfied then 
Equation © or © is used to determine the inward migra- 
tion rate. The rate at which the planet grows is then given 
by Equation (0) with B either a constant (in the "planet 
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dominated" regime) or a constant reduced by 1/B (in the 
"disc dominated" regime). We stop either when the planet 
reaches a — 0.1 AU (the inner edge of our disc models) or 
when the gas disc has dispersed. To synthesise a population 
of exoplanets, we repeat the above for different disc models 
and for different initial starting times and initial semimajor 
axes. 



3 RESULTS 

In these simulations we assume that the sudden change in 
opacity at the snowline (|Menou fe Goodmarfeooll ') results in 
planetary cores being preferentially found beyond the snow- 
line, rather than inside the snowline. We therefore assume, 
somewhat arbitrarily, that the exoplanet density is about 6 
times lower inside the snowline than it is just outside the 
snowline. We also assume that type I migration results in 
a pile-up at the snowline and hence that the distribution 
beyond the snowline is uniform in log a. The only other free 
parameter we have is /?, the ratio of the gas accretion rate 
onto the planet to the gas accretion rate through the disc. 

Figure [4] shows results for 3 different planet growth 
rates, j3. The left-hand panel is for j3 — 0.5, the middle panel 
is for P — 0.7, and the right-hand panel is for P = 0.9. The 
top row of figures shows planet mass against semimajor axis 
for ~ 1000 simulated planets and with the semimajor axis 
normalised with respect to the snowline of the planet's host 
star. For each planet, we have also randomly selected an in- 
clination angle such as to produce an isotropic distribution 
of orbit orientations. We also only show those that would 
have been detected with a 20 year radial velocity campaign 
with a cadence of 1 month and a rms velocity sensitivity of 
lms. The diagonal lines are the same as that in Figure 
[2] Figure [4] already shows that for P = 0.5 the density jump 
occurs inside the diagonal lines, while for f3 = 0.9, it appears 
to be - largely - beyond the diagonal lines. The best match 
to that observed (Figure [2} appears to be the middle figure 
in the top panel in Figure 2] in which f3 — 0.7. 

As discussed in relation to Figure [3] we can use the di- 
agonal lines to predict where each planet emerged prior to 
the start of type II migration and associated gas accretion. 
This is shown in the middle row of Figure [4] The middle 
figure of this row (P = 0.7) is very similar to that seen in 
Figure [3] and suggests that, prior to type II migration, the 
planet population peaks at the snowline and that planets 
then grow in a manner consistent with the diagonal bound- 
ary in Figure[5] However, for j3 = 0.5 and /3 = 0.9, the peaks 
fall inside and beyond the snowline respectively. 

The slope of the diagonal boundary is largely deter- 
mined by the relationship between the rate at which the 
planet migrates inwards and the rate at which it grows. The 
bottom row of Figure [4] shows a sample of growth tracks for 
each /3 value. This illustrates that the slope of the growth 
tracks in mass-semimajor axis space depends on the value of 
P (the ratio of the accretion rate onto the planet to the gas 
accretion rate through the disc). The slope of the boundary 
in mass-semimajor axis space (top figures in Figure ap- 
pears not to precisely match the slope of the growth tracks, 
but is clearly influenced by the slope of these growth tracks. 
If P > 0.7, the diagonal line would be steeper than that 
observed in Figure 1, while if P < 0.7, the diagonal line 



would be shallower. The tracks also show that these planets 
are typically migrating in the "planet dominated" regime 
rather than in the "disc dominated" regime. 



4 DISCUSSION AND CONCLUSIONS 

We use self-consistent disc simulations together with mod- 
els of gas giant planet migration and growth to show that 
if the snowline influ e nces core migration, as suggested by 
Menou fc Goodman] (|2004l ), we can largely reproduce the 
observed planet mass against semimajor axis distribution. 
This in itself is interesting as it indicates that the snow- 
line plays a crucial role in preventing planetary cores from 
migrating (via type I migration) into the host star. Further- 
more, we use this to quantify the gas accretion rate onto the 
planet during the planet dominated type II migration phase. 
We find that it must accrete at a rate of about 70 % that of 
the gas accretion rate through the disc. This is the first time 
that we have been able to quantify the rate at which a gas 
giant planet grows during the final stages of disc evolution. 

Our results indicate that during type II migration a gas 
giant planet consumes ~ 70 % and allows ~ 30 % of the gas 
to flow inward. As a result, we expect a slightly lower gas 
density inside the planets orbit. It has also been suggested, 
howe ver, that dust filtration at the gap edge l|Rice et al.l 
2006) will prevent all but the smallest dust grains from 
reaching the inner disc, significantly enhancing the gas-to- 
dust r atio. This is consisten t with observations of transition 
discs |Espaillat et al. 1120101 ) which have near-IR deficits in- 
dicating a lack of warm dust in the inner disc, but still ap- 
pear to be accreting at TTauri-like rates. 
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Figure 4. Series of figures for 3 different planet accretion rates. The planetary accretion is represented by /3, the ratio of the accretion 
rate onto the planet to the gas accretion rate through the disc. The left-hand panel is for = 0.5, the middle panel is for /3 = 0.7, and the 
right-hand panel is for f} = 0.9. The top row shows the final state distribution of planet mass and orbit semimajor axis, normalised with 
respect to the snowline, after the planets upward growth and inward migration are halted by photo-evaporation of the disc. The middle 
panel is the distribution of the initial semimajor axis (a a determined in the same way as for Figure [3] and normalised with respect to 
the snowline of the host star) for the 3 different boundaries illustrated by the solid, dashed and dotted lines in the top set of figures. The 
bottom is a sample of 100 planet growth tracks showing how the planets grow while migrating inwards through the disc. 
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APPENDIX A: DATA ANALYSIS 

To measure the slope of the observed linear feature in 
log M/Mj-log a/a snow space, we aim to fit a straight line 

log(a/a snow ) =c + a[log(Msini/Mj)-log7], (Al) 
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Figure Al. Plot of Msini against a/a snow for the observed 
planet distribution. The box shows the selection of planets that 
are used to fit a linear density feature. The diagonal line (equa- 
tion [Al} is the best fitting line that maximizes the density ratio 
on each side of the line. Other details are explained in the text. 



to the data so that it maximizes the density ratio on either 
side of the line. This is related to the fitting parameters 
(5,7]) that we use in the paper above through 77 = 1/a and 
S — io c_alos7 . To exclude other features of the simulated 
and observed planet distributions we select only planets with 
0.2 ^ M/Mj < 7 and 0.07 ^ a/a sllow < 1, and measure the 
planet density on each side of this box, as divided by the 
line. In order to penalize solutions that maximized the den- 
sity ratio by producing a small area on either side of the 
line, we actually maximized the quantity ^N+N-p+/p-, 
which we call the fitting metric, where is the planet 

density to the right/left of the line, and N + /_ is the num- 
ber of planets to the right/left of the line. For the observed 
planets, to account for the varying detection efficiency as a 
function of the radial velocity semiamplitude K, we consid- 
ered planets with K < 2.84 m s _1 to contribute a count of 
(2.84 m s^/K) to p +/ _ and N +/ _. 

Figure IA1I shows the observed planet distribution with 
the best fitting line, which has parameters c = —1.005 and 
q = —0.350. It is clear by eye that the number density 
of planets is significantly higher to the right of this line - 
the density contrast is p+/p_ = 3.83. In the plot, plan- 
ets are shaded by their RV semiamplitude, such that the 
planets with the lightest grey are unweighted (i.e. K > 
2.84 m s _1 and they count as 1 planet) and black indicates 
K = m s" 1 . The point with the smallest RV semiampli- 
tude in the fitting box has K — 0.91 m s~ and so counts 
as 3.1 planets. 

Figure IA2I shows the fitting metric ^N+N-p+/p- 
mapped against the two parameters of the line. However, it 
is not possible to estimate the parameter uncertainties from 
the fitting metric as it does not obey \ 2 statistics. Instead, 
we can estimate the uncertainties by a form of bootstrap- 
ping. To estimate the uncertainties by bootstrapping, we 
produced fake planet distributions by Poisson sampling from 
the observed planet number density (including K weighting) 
in four uneven quadrants of the fitting box. The interior 
boundaries of the quadrants were defined by the best fitting 
line and the horizontal dashed line in Figure |AH at the center 
of the selected mass range. The density in each quadrant was 




Figure A2. Map of the fitting metric plotted against the param- 
eters c (intercept) and a (slope) of the line fitted to the density 
feature seen in the observed planet distribution. The cross marks 
the best-fit solution for the observed distribution, and numbers 
denoting the j3 value for simulated planet distributions are placed 
at the corresponding best-fit parameters. The black contours show 
estimated 68- and 95-percent confidence intervals for the observed 
planet solution based on a form of bootstrap resampling, while 
the dashed line was used to discard outlier solutions from this 
bootstrap analysis. 



assumed to be uniform in log(a/o sn ow) and log (A/ sini/Mj). 
Each fake distribution is then fitted in the same way as the 
observed distribution, yielding an estimate of the parame- 
ters c and a in each case. A small number (6 percent) of 
the bootstrap samplings produce solutions far from the best 
fit parameters for the observed planets, and are distributed 
roughly uniformly over the lower right half of the c-a param- 
eter space (in the lower, green region in Figure [A"2)l , while the 
majority of samplings lie near the observed planet solution, 
with a distribution that roughly follows the degeneracy sug- 
gested by the map of the fitting metric. The outlier points 
significantly skew estimates of the error ellipse so we discard 
samples that fall below the dashed line in Figure IA2I The 
contours in the figure show the 49-, 70- and 85-percent con- 
fidence intervals for a Gaussian with mean and covariances 
matching the remaining bootstrap sample, while the cross 
marks the position of the solution for the observed planet 
distribution and numbers indicate the solutions for the sim- 
ulated data sets with different values of j3. It is clear that 
of the simulated planet distributions, the properties of the 
linear density feature in the simulation with j3 — 0.7 best 
match those of the feature in the observed distribution. 
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